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Abstract 

Our previous works have shown the statistical mechanics of self-gravitating 
system. In this paper, we will show its thermodynamics and compare our 
results with observations and simulations. We propose that our statistical 
mechanics can be based on ergodicity breaking and Boltzmann entropy, and 
its assumptions do not contradict with the reality. With the principle of sta- 
tistical mechanics, wc will show our definition of temperature and then the 
capacity of self-gravitating systems. We find that the gravothermal catas- 
trophe may be a special case of our theory. Our results also provide new 
explanations for the density profiles of observations and numerical simula- 
tions, especially we think that the non-universal density distribution in the 
simulations of dissipationless collapse is not caused by the different initial 

collapse factor. We will also discuss about the core-cusp problem. 
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1. Introduction 

Statistical mechanics of systems with long-range interactions presently 
has not been generally considered to be well established, because thermody- 
namics of long-range system is very different: the energy is non-additive; the 
entropy needs not to be a concave function of energy^ so the specific heat may 



2|, |3| ; then the microcanoni- 



be negative in the microcanonical ensembles 
cal ensembles and canonical ensembles may be not equivalent; the ergodicity 
which is the fundamental assumption of statistical mechanics may be broken 
; and many others. 

However, our latest studies provide different understandings of the long- 



range statistical mechanics. 



5|] preliminarily studies the fluid entropy and 



proposes a self-gravitating system's entropy principle which is different from 
the well known principle of maximum entropy. 6|] completes the variation 
process of entropy and confirms that this entropy solves the problem that 
people always obtain the infinite mass and energy when using the Boltzmann- 
Gibbs entropy. Now it is necessary to show the thermodynamics correspond- 
ing to our statistical mechanics. 

Besides, presently there have been many observations and increasing nu- 
merical simulations of some astronomical systems which can be approxi- 
mately considered to be self-gravitating. While there are many observations 
of globular clusters and elliptical galaxies, there are few simulations of dis- 
sipationless collapse of galaxies which will interest us here. More current 
simulations are cosmological simulations which have been successful to be 



3 



consistent with the universe's large scale structure predicted by the ACDM 
model, but there seems to be some contradictions with observations at the 
galaxy scale, such as that current cosmological simulations provide us a NFW 



density profile 



i 



which has a central cusp and is not consistent with the measurements of Low 
Surface Brightness galaxies jsl that show us a central core. Other two incon- 
sistencies between simulations and observations include the missing satelhte 



and the problem about the distribution of angular momentum 26|]. In this 
letter we will also discuss some of these results. 

2. Equation of state 

We first summarize some of our previous works with more observational 
evidences. Our attentions are restricted to the spherical self-gravitating sys- 
tem. It is known that if we insist on the ergodicity and the principle of 



maximum prmcip 
mass and energy 



we will get the isothermal solution which gives infinite 



10|. While Mukamel ^ has noticed that the ergodicity is 



easily broken for long-range systems, which means that not all the microstates 
that satisfy the macroscopic constraints can occur in the long-range systems. 
From observations and simulations flQ we know that the velocity disper- 
sion (the "temperature") is different at different places of the system, while 
the velocity distribution at each point is nearly Gaussian [isj]; besides, the 
force on the particle is mainly determined by the large scale structure of the 



4 



self-gravitating system lOj , which may indicate that the particles in a micro 
volume d^x can be treated as being free. Based on above conditions, we 
assume that the phase space distribution can be written as 

/ oc e (2) 



where = P/ p and P is the effective pressure defined in 



^ = ^ + 2/3i^. (3) 
dr dr r 

Pi is the parameter of velocity anisotropy, and pr = P^r- We think that 
eq. ([2]) may mean that the ergodicity is broken for the whole system but may 
still set up locally for arbitrary micro volume d^x, because if the system's 
ergodicity can set up, we can derive the Maxwellian velocity distribution 
from statistical mechanics, while eq.([2]) means that the Maxwellian velocity 
distribution only can set up in the micro volume d^x. Then from the form 

n 

of Boltzmann entropy we can obtain (also see section 4.3 of |jM] for the 
calculations) : 

S= //ln/dr = 47r- / pln(^)rMr, (4) 

m J m Jo p°' 

where p is the space density, k is the Boltzmann constant, and m is the mass 
of single particle. Notice that jisl also studies the elliptical galaxy (which 
is commonly treated as being coUisionless) by using the ideal gas's entropy 
which is almost the same as eq.Q. We think that eq.(l2]) is applicable both 
for the systems because eq.(|2]) depends only on (j^ which can exist even that 
the two-body relaxation can be neglected or not; and whether the two-body 



relaxation can be neglected or not, when the system evolves slowly and is 
close to static, we can both obtain Q 

— = -P— (5) 
dr dr 

But the difference is that if the two-body relaxation can not be neglected. 



the / in the eq. 



can be fine-grained; while for the coUisionless system, the 



Jeans theorem [10| does not allow the form of eq. ([2]) , so we need to let / to be 



coarse-grained, and another reason to do this is that S can not evolve if S is 

n 

a functional of the fine-grained phase density [17J. Then with the constraints 
of mass and energy 

/•oo 

M = 47rpr^dr, 
Jo 

/•oo roo 

E = —2710 I p{r)m{r)rdr = —671 / Pr^dr, (6) 
Jo Jo 

and the differential constraint provided by eq.(l5]), we can calculate the en- 
tropy's extremum which requires 

SSt = 5(3 /k) - I35E - a6M -5 [ rfrr/(r)(^ + p^) = 0, (7) 



dr dr 

where a, /3 and T]{r) are Lagrangian multipliers. Notice that because of 
the definition of P we do not require that the velocity anisotropy to be 
zero during the process of calculating the extremum of the entropy, while 
observations and simulations also show the non-zero anisotropy in the outer 
region. Eq.(I7])'s solution can be approximately written as an equation of 
state jsl 

p = m(3P + aP"", (8) 
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Where n=3/5 or 4/5. The second order variation of St is 

52S,^i/[(4vr + l)/3^)(5M)2 
M dr, 

which indicates that St is a saddle point (/3 > 0). This agrees with the con- 



clusion of 



3. 



We can understand it by that based on the broken ergodicity 
the equilibrium state of self-gravitating system corresponds to the time when 
the number of microstates is an extremum and does not need to be a maxi- 
mum , so S" is allowed to be not a maximum. Notice that the existence of a 
ensures that the system can have finite mass and finite energy. High a means 
small extent (radius) of the system and High /3 requires a small central core 
(see Fig.E]). 

3. Thermodynamics of self-gravitating systems 

In the classical thermodynamics of self-gravitating systems, because of the 
virial theorem, 2K + W = where K and W are the kinetic and potential 
energy respectively, E = K + W = —K. While people define the system's 
temperature by = 3NkT/2 where N = M/m, so 

E = --NkT, = _ < 0. (10) 

An important conclusion is that the negative capacity can cause the gravother- 



mal catastrophe because of the principle of maximum entropy [10|]. However, 
here we argue about the definition of the temperature K = 3NkT/2 which 
has not been proven to be applicable to self-gravitating systems: from eq. fllOl) 
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Figure 1: The density profiles with different values of {3 and fixed a. po is the central 
density, is a characteristic scale, and we assume that these two quantities have been 
given. With the increase of (3 from to oo, the region of the central core becomes small 
and the density profile tends to the singular isothermal sphere. 
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we find that E oc N for fixed T, which is not consistent with the non-extensity 



of long-range systems [19[; although we show that the heat capacity is neg- 
ative, but it is still positive when the density contrast between the center 



and the outer part is small enough 



10|. Based on above, we may need to 



reconsider the thermodynamics of self-gravitating system. 

According to the method of statistical mechanics, in the eq.(I7]) we can 
identify 

/3 = 1/kT. (11) 

This means that we have found a self-gravitating system's thermodynamic 
equilibrium state characterized by the temperature T, which is a very differ- 
ent conclusion from all previous works. Then a natural question is how to 
reach this equilibrium state, which will require us to consider the capacity. In 
the following, we first take n = 3/5. From figure. [2], we know that commonly 
eq.® can be approximated as kTp/m = P, which will be substituted into 
the second term of eq. ([H]) 

kT , , .33. . , 

p=_ _ 5p5 , 12 

m m 

Which is reminiscent of the van der Waals equation that is to describe the 
non-ideal gas with the interactions among molecules: 

2 

Ti n 

{p + ^){V-nb) = NkT, (13) 

where a and b are related to the attractive and repulsive forces respectively. 
If we take p = Nm/V and set 6 = 0, we find that eq.f ll3p will be very similar 
to eq.( fT2|) . which may further support the rightness of eq.( fT2|) . Multiplying 



eq. (fT^ with —3/2 and then integrating the both sides of the equation over 
the whole volume, we can get 

E = --kTlN -a{kT)lA],A = An H r^-^f/^dr. (14) 
2 Jo m 

Notice that here E is not proportional to A^, which just is the requirement 
of long-range systems. Bounded systems satisfy E < 0, which requires that 

a{kT)l < ^. (15) 

Then we find that the heat capacity is 

Cv=g;^ = ^k[-a{kT)lA-N], (16) 

which can be positive or negative. However, according to our results the 
self-gravitating system can always approach the thermodynamical equilib- 
rium whether the capacity is positive or not, which will be analyzed as the 
following: 

(1). when C > 0, we can easily obtain 

'4 < ami < ^. (17) 

8A ^ ^ A ^ ^ 

For a given value of a, C > is caused by high temperature, then from eq.(j9]) 
we know that the (5M)^ term can be neglected, so the entropy is maximized 
and the thermodynamical equilibrium can be approached. This also can be 
understood by that when T is high the self-gravity of our system is very week 
so that we can treat the system as an ideal gas. When T tends to be infinite, 
the density profile is a homogeneous distribution truncated at a finite radius, 

10 



just as shown in figured This truncation is caused by the existence of the 
Q,p3/5 f^QYm, which shows us that the existence of a just is equivalent to add 
an insulting shield to protect the particles from escaping and to ensure of 
the conservation of the total energy of the system. 
(2). C < requires that 

< a(A;T)3 < — . (18) 

In this case, the classical thermodynamics tells us that the system will be- 
come instable and result in gravothermal catastrophe. However, in our new 
thermodynamics, the entropy is a saddle point and is not maximized, 
which means that the second law of thermodynamics is not valid for self- 
gravitating systems. We first show an extreme case: when T — > 0, from 
eq.([n]) we know that the (5M)^ term becomes so large that S is minimized, 
so considering the Clausius statement of the second law of thermodynamics, 
the heat now is allowed to spontaneously be conducted from low temperature 
system to high temperature system (In logics we know that if the original 
proposition is correct, then its converse-negative proposition will also be cor- 
rect. While current original proposition can be that "if the second law of 
thermodynamics is correct, then the entropy of the equilibrium state will be 
maximized" , its converse-negative proposition is "if the entropy of the equi- 
librium state is not maximized, then the second law can be not correct, i.e. 
the heat now is allowed to spontaneously be conducted from low tempera- 
ture to high temperature"), and the thermodynamical equilibrium can still 
be approached although C < 0. We also notice that in this case the density 
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profile approaches the singular isothermal sphere and self-gravity becomes 
very strong, so the self-gravity is the origin of the negative heat capacity, 
which is consistent with previous study [2,]. In the common cases the entropy 
is a saddle point, and the thermodynamical equilibrium can be approached 
based on above analysis. 

Before comparing our results with observations and simulations, we first 
talk about the value of n. In 6|] we have said that n = 3/5 in fact is to treat 
$ = for large enough r, so $ = —GM/r + GM/R for r < R where R can 
be the radius of the system, and $ will be zero for r > R; while n = 4/5 is 
to treat $ = —GM/r even for very large r. Evidently for observations and 
simulations the latter one is more accurate than the former one, which also 
can be seen in Figj2l 

4. Discussions for results of observations and simulations 

In 6|] we have showed that eq.® is very consistent with observations. 
It can be approximated as truncated cored isothermal sphere which just 
is the King model I^Qj] that describes the density distributions of globular 
clusters. Besides, although the surface brightness profiles of elliptical galaxies 
are described by the i?^/^ law, some standard elliptical galaxies, such as NGC 
3379, still have a central core, which makes eq.® more consistent with their 



surface brightness than i?^/^ law 
However, some observations 



2l| and simulations 22| have shown that 
many globular clusters are "core collapsed" which commonly is treated as 
an evidence of gravothermal catastrophe but can be stopped by the binary 
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stars whether they are primordial or produced by gravothermal catastrophe 



10|. But our understanding is that because of not vahdity of principle of 
maximum entropy the gravothermal catastrophe may not exist: the den- 
sity profile of core collapsed clusters 10| tends to be the singular isothermal 



sphere, so according to the analysis of above section the initial conditions 
(mass,energy,etc) of these clusters make them have strong self-gravity so 
that their final temperature T — > (notice that the these clusters are closer 
to the galaxy center, which means that they can survive under the strong 
tidal force of the galaxy, so it is reasonable to believe that they have strong 
self-gravity), and they are just approaching or have reached the thermody- 
namical equilibrium. Evidently our explanation does not contradict with 
the observations and simulations that may manifest the gravothermal catas- 
trophe, but there is another result that may more support our theory: the 
gravothermal catastrophe will occur once the density contrast between the 



2|], while 2l| finds that there 



outer and the inner parts is larger than 708.61 
are a number of globular clusters that have very dense core and very short 
relaxation times, but they seem not suffering core collapse and are well fitted 
by the King models. 

Then we should compare our results with the simulations with isolated 



boundary conditions. [23| simulates the dissipationless collapse of galaxies 
and their figure. 4 shows that for the initial uniform spheres the final density 
profiles do not have a universal shape, which is thought to be caused by the 



galaxies' initial collapse factor 2K/W. But in our view, the results of 



23| 



just perfectly confirm the rightness of eq.([H]): its all kinds of density profiles 
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Figure 2: van Albada 
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logr 

's results fitted by eq.®. Q = 2K/W. 



still can be exactly fitted by eq.dH]) with different choices of (3, a and p(0) 
especially when n = 4/5, which is shown in Fig|2l So we conclude that the 
non-university of the final space density may be not caused by different initial 
collapse factor but is caused by different total masses and total energies of 
the simulated systems. 

Until now, our results are consistent with observations of star clusters 
and elliptical galaxies and simulations of dissipationless collapse. Evidently 
our results can not explain the cusp profile which is to describe the density 



profile of the simulated dark matter halos. But we also notice that: 2^ also 
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Figure 3: Density profiles of dissipationless simulations with different initial conditions, 
compared with eq.® and NFW profile. R^q is the radius with half of mass. We think 
that the deviations from NFW profile at the outer part is caused by that CDM halos are 
not strictly isolated. 
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simulates the above collapse with many kinds of initial conditions (but all 
with the isolated boundary condition): their figure. 3-8 show us that for 
most initial conditions, the final density profiles have a central core and also 
can be exactly described by eq.([H]); but for the clumpy initial distribution, 
the cusp can be easily produced and the final central density is better fitted 
by eq.([T]). All above have been shown in FigJS] Besides, if the initial den- 
sity slope is very steep, a cusp will also be produced, but not as quickly as 



25| that performs the cosmological 



the clumpy case. While Syer & White 
simulation has ever suggested that the cusp profile is a result of hierarchy of 
clustering by mergers of the smaller clumps to the bigger clumps. However, 



Moore et al [26|] has almost rejected this possibility by the simulations of 
warm dark matter which can smooth the density perturbation at the galaxy 
scale. So we can conclude that the cored profile only can appear in the simu- 
lations with isolated boundary condition and homogeneous initial condition 
at the galaxy scale, and any disagreements with one of these two conditions 
will lead to the cusp profile. Besides, observations of LSB galaxies and others 
strongly support the cored profile; the similarity solution of the secondary 
infall model of hierarchy clustering, which predicts that the final density pro- 
file depends on the initial density perturbation, has shown us that for the 
initial density perturbation 5i ~ M% the resulting density slope is 

dlnp 9e , , 



which is a cusp 27| ; our results about the thermodynamical equilibrium state 



of an isolated self-gravitating system suggest the cored one. If we neglect 
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the interactions between baryons and dark matter halos, do these evidences 
indicate that the galaxies formed with more isolated boundary condition and 
more homogeneous initial condition? 

5. Conclusion 

In conclusion, if we assume that the final state of self-gravitating system 
is the state when the number of microstates (so the entropy) is an extremum, 
which does not contradict with the ergodicity breaking, we can find a different 
thermodynamics of self-gravitating system. The equation of state of the 
equilibrium state is similar to the van der waals equation, and the final 
density profile is determined by M and E which are used to control the 
value of P and a. 

Our different thermodynamics states that the equilibrium may always be 
approached. When the temperature is large enough, our thermodynamics 
can come back to the classical thermodynamics, and the equilibria of self- 
gravitating system can be described by the isothermal sphere; when the 
temperature tends to be zero, the heat capacity is negative, but the the 
entropy also is minimized, so the system may still can be in equilibrium and 
gravothermal catastrophe may be not necessary to exist. Some observations 
and simulations that manifest the gravothermal catastrophe can be treated 
as a special case of T — > in our results. 

When we pay attention to the simulations of dissipationless collapse, we 
find our results can exactly fit their results and show that the the non- 
universal final density distribution may be not caused by the collapse factor. 
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We also discussed about the core-cusp problem, which may cause many ef- 
fects at current understandings of galaxy formation. Our results and some 
explanations are really different from previous works, which needs to be fur- 
ther confirmed. 
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